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At s tract 


The dynamical behavior of spacecraft propellant affected by the asymmetric 
combined gravity gradient and jitter accelerations, in particular the effect of 
surface tension on partially- filled rotating fluids applicable to a full-scale 
Gravity Probe-B Spacecraft dewar tank has been investigated. Three different 
cases of orbital accelerations : (a) gravity gradient-dominated, (b) equally 
weighted between gravity gradient and jitter, and (c) gravity j itter- dominated 
accelerations are studied. The results of slosh wave excitation along the 
liquid-vapor interface induced by gravity gradient- dominated accelerations 
provide a torsional moment with tidal motion of bubble oscillations in the 
rotating dewar. The results are clearly seen from the twisting shape of the 
bubble oscillations driven by gravity gradient-dominated acceleration. The 
results of slosh wave excitation along the liquid-vapor interface induced by 
gravity jitter-dominated acceleration indicate the results of bubble motion in 
a manner of down-and-up and leftward- and- rightward movement of oscillation when 
the bubble is rotating with respect to rotating dewar axis. Fluctuations of 
angular momentum, fluid moment and bubble mass center caused by slosh wave 
excitations driven by gravity gradient acceleration or gravity jitter 
acceleration are also investigated. 
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I. Introduction 


In space experiments, superconducting sensors are used for gyro read-out 
and maintain a very low temperature for mechanical stability. The approaches to 
both cooling and control involve the use of superfluid liquid helium. For 
example, to test Einstein's General Relativity, the Gravity Probe-B (GP-B) 
Spacecraft adopts the boil-off from the cryogenic liquid helium dewar as a 
propellant to maintain the attitude control and drag- free operation of the 
spacecraft 1 . The potential problems for cryogenic liquid in the dewar container 
could be due to asymmetry in the static liquid helium distribution which is 
controlled by the capillary effect governed liquid-vapor interface disturbances 
in a microgravity environment. 

In the absence of temperature gradient along the liquid-vapor interface for 
superfluid helium which drives Marangoni convection, the equilibrium shape of the 
interface in a rotating container is governed by a balance of capillary, 
centrifugal and gravitational forces in microgravity. The determination of 
liquid-vapor interface profiles based on computational experiments can uncover 
details of the capillary effect driven flow which can not be easily visualized 
or measured experimentally in a raicrogravity environment. 

The instability of the capillary effect governed liquid-vapor interface 
profiles can be induced by the presence of longitudinal and lateral 
accelerations, vehicle vibration, and rotational fields of the spacecraft. Thus, 
slosh waves are excited, producing high and low frequency oscillations in the 
fluid systems. The sources of the residual accelerations range from the effects 
of the Earth's gravity gradient acceleration and jitter acceleration, which 
includes atmospheric drag on the spacecraft, spacecraft attitude motions arising 
from machinery vibrations, thruster firings, crew motion, etc. A recent study 2 
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suggests that the high frequency accelerations may be unimportant in comparison 
to the residual motions caused by low frequency accelerations. 

The time -dependent dynamical behavior of the capillary effect governed 
liquid-vapor-solid interface disturbances for rotating fluids in reduced gravity 
environments was simulated by numerically solving the Navier Stokes equations 
subject to the initial and boundary conditions 3 ' 4 . At the interface between the 
liquid and the vapor fluids, both the kinematic surface boundary condition and 
the interface stress conditions for components tangential and normal to the 
interface were applied 3 " 6 . The initial conditions were adopted from the 
steady- state formulations developed by Hung et al 3 with an example in the 
geometry of the GP-B Spacecraft 1 . Some of the steady-state formulations of 
interface shapes were compared with the available experiments carried out by 
Leslie 7 in a free-falling aircraft (KC-135) . The experiments carried out by 
Mason et al s showed that the classical fluid mechanics theory is applicable for 
cryogenic liquid helium in large containers. 

As for the effect of centrifugal force on the liquid helium and helium 
vapor with an example of the GP-B Spacecraft, a spin rate of up to about 1 rpm 
will be imposed for instrument calibration in the early stages of the experiment. 
After calibration, the rotation rate will be reduced to its operational value of 
approximately 0.1 rpm. 

As the spacecraft moves along the orbit, any fluid mass capable of motion 
relative to the spacecraft is subject to the acceleration that arises from the 
gravity gradients of the Earth 9 11 and orbital environmental gravity jitter 
acceleration. In this paper, in addition to the ubiquitously available study of 
sloshing dynamics modulated liquid-vapor interface oscillations in normal and 
microgravity environment, we are particularly interested in the investigation of 
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orbital acceleration, the combined gravity gradient and jitter accelerations, 
induced sloshing dynamics and the characteristics of its impact on orbital 
spacecraft fluid systems. In this study, comparisons are made of the dynamical 
evolution of sloshing dynamics excited moment and angular momentum fluctuations 
of fluid system for the combined gravity gradient and jitter accelerations at 
different background of gravity jitter accelerations due to the various ranges 
of environmental drag. In other words, the study is based on the orbital gravity 
gradient acceleration with a combination of three ranges of gravity j itter 
accelerations (10 6 , 10 7 and 10~ 8 g 0 ) for the purpose of studying the dynamics 
of capillary effect governed liquid-vapor- solid interface disturbances and their 
effects on fluid moment and angular momentum disturbances. 

II. Mathematical Model of Capillary Effect Governed Slosh Wave Excitation 

Along the Liquid-Vapor Interface 

Consider a closed circular dewar cylinder of an inner radius a lf and an 
outer radius a 2 , with height L, which is partially filled with cryogenic liquid 
helium, and rest of the ullage is filled with a helium vapor. Density and 
viscosity of liquid helium and helium vapor are p L , /j l , p v , and respectively. 
Let us use cylindrical coordinates (r, 6 , z), with corresponding velocity 
components (u, v, w) , and corresponding residual gravity acceleration, such as 
gravity gradient components (a ggr , a gg<? , a gg 2 ) and gravity jitter components 
( a sj,r> a gj,0> a gj,z)* The whole fluid system is spinning in the z-axis with an 
angular velocity of u> and an angular acceleration of a. The governing equations 
for non- inertia frame bound coordinates are illustrated in our recent studies 12 ' 13 
and will not be repeated in this paper. In these formulations, 2pu>u and 2ptov 
denote Coriolis forces; pro, the angular acceleration; and pro; 2 , the centrifugal 
force of spinning motion. 
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In order to solve sloshing dynamics problems of liquid propellant systems 
in orbital spacecraft under a microgravity environment, one has to solve the 
governing non- inertia frame equations 12 ' 13 accompanied by a set of initial and 
boundary conditions. A detailed illustration of these initial and boundary 
conditions concerning the sloshing dynamics of fluid systems in microgravity were 
precisely given in our recent studies 13 ' 14 and will not be repeated in this paper. 
In this study, in order to show a realistic example, a full scale GP-B spacecraft 
propellant dewar tank with an inner radius of 12 cm and an outer radius of 68 cm 
and a height of 145 cm has been used in the numerical simulation. The propellant 
tank is 80% filled with cryogenic liquid helium and the rest of ullage is filled 
with helium vapor. The temperature of cryogenic helium is 1.8 K. In this study 
the following data was used: liquid helium density - 146 kg/m 3 , helium vapor 
density - 1.47 kg/m 3 , fluid pressure - 1.66 Pa, surface tension coefficient at 
the interface between liquid helium and helium vapor - 0.0353 N/m, liquid helium 
viscosity coefficient - 9.61 x 10~ 9 m 2 /s ; and contact angle - 5°. The initial 
profiles of the liquid-vapor interface used in this study have been given 
explicitly through the steady state computations made by Hung and Leslie 15 and 
Hung et al. 3 ' 4,16 , which were checked by experiments carried out by Leslie 7 . 

A staggered grid for the velocity components is used in this computer 
program. The marker-and-cell method of studying fluid flows along a free surface 
is adopted 17 20 . The formulation for this method is valid for any arbitrary 
interface location between the grid points and is not limited to middle point 
interfaces 19 . An algorithm for a semi - implicit method 20 was used as the procedure 
for modeling the flow field. The time step is determined automatically based on 
the size of the grid points and the velocity of flow fields. More than 200 times 
of iteration have been applied for each time step to assure the convergence of 
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series. As for the volume conservation of the liquid, a deviation of less than 
1% error of volume is always guaranteed before a move to the next time step. 

III. Dynamical Behavior of Liquid-Vapor Interface Disturbances Driven by 
Orbital Accelerations in Microgravity 
For the case of the GP-B Spacecraft, which Is an Earth satellite orbiting 
at 650 km altitude directly over the poles, the orbit period, r Q can be computed 
from the following expression: 


x 0 =2k 


R 


3/2 

c 




( 1 ) 


where R E denotes radius of Earth (- 6373 km); R c , radius of circular orbit (- R E 
+ h — 7023 km); h, orbit altitude (— 650 km); and g Q , Earth gravity acceleration 
(*• 9.81 m /s 2 ). For the case of GP-B, the orbit period r 0 — 97.6 min, and orbit 
rate n - 2n /r 0 - 1.07 x 10~ 3 rad / s . 

The gravity gradient acceleration acting on the fluid mass of spacecraft 
can be shown as 


S gg = -n 2 [3 (f c ‘d) f c -d] (2) 

where a gg denotes gravity gradient acceleration vector; 3, vector from the fluid 
element to the spacecraft mass center; r c , unit vector from the spacecraft mass 
center to the center of the Earth; and n, the orbit rate. 

For the case of the GP-B Spacecraft, it is assumed that the gravity exerted 
on the mass center of the spacecraft orbiting around the Earth on its specified 
orbit is zero. In other words, all the gravity acceleration exerted on the 
spacecraft is nothing but the gravity gradient acceleration which is defined in 
Equation (2). In this study, we are interested in investigating how combined 
gravity gradient and jitter accelerations affect the dynamical behavior of 
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cryogenic fluid elements of helium. 

For the convenience of mathematical calculation, let us describe all the 
parameters involved in Equation (2) in terms of cartesian coordinates. In order 
to match with the computer simulation, mathematical derivation are considered in 
the first quadrant. Figure 1 illustrates the geometrical relationship of the 
parameters shown in Equation (2) . 

Let us consider the fluid element of interest, m, located at (r, 0, 2 ) in 
cylindrical coordinates and at (x, y, z) in cartesian coordinates. As |d| is 
much smaller than the distance between the location of the GP-B Spacecraft 
geometric center to the center of the Earth, r c through the GP-B geometric center 
and r c through the fluid element, m, is basically the same. 

The gravity gradient acceleration acted on the fluid element located at (r, 
0, z), in non-inertia frame coordinates is shown as follows 12 ' 13 ' 16 : 

a 1 cosO+cot) sin(0+a>t) 0 

S gg = a gg,e = - n 2 -Sin(0 + CDt) COS (0+0) t) 0 
a gg, r 0 0 1 

3 [ (z-L/2) costy E -rcos (0+cot) *simJ/ £ ] sini}r £ +rcos (0+co t) 

rsin(0+cot) ( 3 ) 

-3 [ ( z-L/ 2 ) cosiJj £ -rcos (0+c*> t) *sirujr £ ] cos t £ + (z-L/2) ] 

where denotes the azimuth angle of Earth toward the location of spacecraft 

geometric center, or the angle between the spacecraft rotating axis and the 

radial direction from spacecraft geometric center to center of Earth; n, the 

orbit rate; and t, the time of spacecraft orbiting around the Earth. 

At time t - 0, the rotating axis of the spacecraft is aligned with the 
radial direction of Earth center to the spacecraft geometric center. An example 
is given by assuming that the spacecraft axis of rotation is linearly turning 
around 0° to 360° in 1200 s, which is twice the period of spacecraft rotating 


6 



time (0.1 rpm is equivalent to a 600 s rotating period), when the spacecraft is 
orbiting around the Earth. Thus, ^ can be defined as 

♦*=2-2.t (4) 

where r is the spacecraft gravity turn-around time (- 1200 s in this case); and 
t is the time measured from the instant when the direction of the spacecraft 
rotating axis is aligned with the radial direction of the spacecraft mass center 
to the center of the Earth. 

The fluctuations of residual gravity due to orbital environmental gravity 
jitter acceleration are modeled by the following equation^: 

^=9 r fl tl + -|- sin ( 27tft ) 1 ( 5 ) 

where g B denotes the orbital background gravity environment, and f(Hz) stands for 
the frequency of gravity- j itters . 

In this study, three ranges of gravity jitter background with gravity of 
10 6 , 10 7 and 10 8 g Q due to different intensity of spacecraft atmospheric drag 
and rotating speed of 0.1 rpm for normal GP-B spacecraft operation, and three 
frequencies of gravity jitter with the values of 0.1, 1.0 and 10 Hz have been 
considered for the investigation of the oscillations of the liquid-vapor 
interface. The components of gravity jitter acceleration in the non- inertia 
coordinate system shall be given by 

£gj = ( a gj , i > a gj , q ' a gj.z"> = [ gs iruj/^os (0+g) t) , in ( 0 + 0 ) C) , -gcosilfjj] 

( 6 ) 

The characteristics of gravity gradient and jitter accelerations, shown in 
Equations (3) and (6), respectively, are quite different. Gravity gradient 
acceleration shows the following characteristics: (a) acceleration acts on any 
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fluid mass inside the container which increases two units of acceleration per 
unit of distance measured from the container mass center (point 0 g at Figure 1) 
to the location of the fluid mass parallel along the radial axis from the 
container mass center to the center of the Earth (parallel to unit vector r c 
shown in Figure 1); (b) the acceleration acts on the fluid mass which decreases 
one unit of acceleration per unit of the shortest distance measured from the 
location of the fluid mass to the radial axis along the vector from the container 
mass center to the center of the Earth 10 . For example, Figure 2 shows the time 
variation of gravity gradient accelerations for a turn-around period of 1200 s 
with a container rotating speed of 0.1 rpm for components along (x, y , z) 
directions acted on the fluid mass located at (r, 6 , z) - (40 cm, jr/4, 10 cm). 
As the magnitude and direction of gravity gradient acceleration acted on each 
fluid mass is strongly dependent upon how far the location of the fluid mass 
deviates from the container mass center measured along the axis parallel to the 
vector r c which varies with respect to time, it can be anticipated that the 
gravity gradient acceleration acted on the fluid mass is different for fluid mass 
at different locations in the container. Figure 2 shows that the magnitude of 
gravity gradient acceleration is on the order of 10~ 7 g Q . On the contrary to 
gravity gradient acceleration which acts different acceleration on the fluid mass 
at different locations in the container, gravity jitter acceleration drives the 
same acceleration on the fluid mass at different locations in the container. 
Figure 3 shows the time variation of gravity jitter accelerations for a turn- 
around period of 1200 s with a container rotating speed of 0.1 rpm and a jitter 
frequency of 0.1 Hz for components along (x, y, z) directions acted on the fluid 
mass everywhere in the container. It is shown in Figure 2 that gravity gradient 
accelerations are on the order of 10‘ 7 g 0 for the spacecraft orbit chosen in 
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this example, three ranges of background gravity, 10‘ 6 , 1CT 7 and 10‘ 8 g,, for 
gravity jitter accelerations correspond to the magnitude of acceleration higher, 
equal, and lower, respectively, than that of the gravity gradient acceleration 
acted on the fluid system of the spacecraft. 

The equilibrium shape of the liquid-vapor interface for a rotating dewar 
under a residual gravity environment below 10'^ and rotating speed of 0.1 rpm 
is very much alike and is a doughnut configuration with a near circular kidney- 
shape cross - section based on the computation of the numerical algorithm developed 
in our earlier studies 3 . Figure 4 shows the initial shape of the interface in 
the r-z plane at 8 — 0° and 180°; Figure 4 shows the initial profile of the 
liquid-vapor interface in the r-z plane at 0 - 90° and 270°; Figure 4 shows the 
initial profile of the liquid-vapor interface in the r-8 plane at height z - 108 
cm; and Figure 4 shows the initial profile of the three-dimensional liquid-vapor 
interface. Because of page limitation, the time evolution of the three- 
dimensional liquid-vapor interface profiles only will be illustrated in this 
paper . 

In this study, the combined gravity gradient and three ranges of jitter 
(10" 6 , 10~ 7 , and 10' 8 g Q ) accelerations acting on the liquid-vapor interface 
oscillations have been investigated. 

( III - A) Combined Gravity Gradient and 10' 8 g 0 Background Jitter Accelerations 
Acting on Liquid-Vapor Interface Oscillations 

In this case, the combined effects of gravity gradient and jitter 
accelerations are completely dominated by the gravity gradient acceleration when 
the 10' 8 g Q and less background gravity jitter accelerations are applied 
simultaneously with gravity gradient acceleration (which is on the order of 10' 7 
go as shown in Figure 2) for the excitation of slosh waves along liquid-vapor 
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interface oscillations. 


Figure 5 shows the time sequence evolution of the liquid-vapor interface 
profiles driven by these combined accelerations. For the convenience of 
comparison, figures of liquid-vapor interface profiles with the same values of 
the time sequences chosen throughout this paper are at time t - 191, 354, 380, 
431, 503, 603, 825, 980, 995, 1050, 1080 and 1200 s. It clearly shows that there 
is a series of asymmetric oscillations excited along the surface of liquid-vapor 
interface driven by asymmetric gravity gradient-dominated acceleration. 

Careful examination of the case of gravity gradient-dominated acceleration- 
effected liquid-vapor interface oscillations , gravity gradient acceleration shown 
in Equation (3) , indicate that there are greater negative components of 
acceleration longitudinal to the direction aligned with spacecraft mass center 
to the Earth center, and smaller positive components of acceleration transverse 
to this direction. As we indicated in Equation (4), angle ife varies with time. 
This phenomena shows that the gravity gradient acceleration exerted on the 
spacecraft is equivalent to the combination of time dependent force with turn- 
around direction and torsional moment acting on the spacecraft when it is 
orbiting around the Earth. It shows in Figures 4 and 5 that the deformation of 
the bubble is created by asymmetric torsional moment with a twisting force. 

Figure 5 shows the time sequence evolution of the three-dimensional 
dynamical behavior of the interface oscillations driven by gravity gradient- 
dominated acceleration. It indicates that the bubble (helium vapor) 
configurations change from axial symmetric to asymmetric profiles at a plane 
aligned with the vector of gravity gradient acceleration. The gravity gradient- 
dominated acceleration produces a combination of time dependent torsional 
movement with tidal motion of bubble oscillations in the rotating dewar when it 
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is orbiting around the Earth. 

(III-B) Combined Gravity Gradient and 1CT 7 g,. Background Jitter Accelerations 
Acting on Liquid-Vapor Interface Oscillations 

In this section, the combined effects of accelerations are equally weighted 
by both the gravity gradient and jitter accelerations when 10’ 7 ^ background 
gravity jitter acceleration is jointly applied simultaneously with gravity 
gradient acceleration (which is on the order of 10~ 7 ^ as shown in Figure 2) for 
the excitation of slosh waves along the liquid-vapor interface oscillations. 
Figure 6 shows the time sequence evolution of the dynamical behavior of the 
liquid-vapor interface oscillations driven by these combined accelerations. 
Similar to Figure 5, the same values of the time sequences are chosen. It also 
shows that there are a series of asymmetric oscillations excited along the 
surface of liquid-vapor interface driven by asymmetric equally weighted gravity 
gradient and jitter accelerations. 

As we recorded in Figure 5, the effect of gravity gradient-dominated 
acceleration exerted on the spacecraft is equivalent to the combination of time- 
dependent force with turn-around direction of torsional moment acting on the 
spacecraft when it is orbiting around the Earth. On the other hand, it will be 
show later in Section (III-C) for the case of gravity jitter-dominated 
acceleration exerted on the spacecraft which is equivalent to time -dependent 
oscillatory forces that push the bubble in the combined directions of down-and-up 
and leftward-and-rightward as the bubble is rotating with respect to the 
spacecraft rotation axis. Figure 6 exactly illustrates these combined effects 
of equally weighted gravity gradient and jitter accelerations which drive the 
bubble in a manner of combined oscillations of the twisting deformation of the 
bubble associated with down-and-up motions. 
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A combined oscillation of twisting bubble oscillations with down-and-up 
bubble oscillations modified by leftward- and- rightward oscillations indicates 
very important characteristics of combined gravity gradient and jitter 
accelerations of equal weight acting on the spacecraft when it is orbiting around 
the Earth. 

(III-C) Combined Gravity Gradient and 10~ 6 ^ Background Jitter Accelerations 
Acting on Liquid-Vapor Interface Oscillations 

In this case, the combined effects of gravity gradient and jitter 
accelerations are completely dominated by the gravity jitter acceleration when 
the 10" 6 g Q and higher background gravity jitter acceleration are applied 
simultaneously with gravity gradient acceleration (which is on the order of 1CT 7 
gQ as shown in Figure 2) for the excitation of slosh wave along the liquid-vapor 
interface oscillations . 

Figure 7 shows the time sequence evolution of the dynamical behavior of the 
liquid-vapor interface oscillations driven by gravity jitter-dominated 
acceleration of 0.1 Hz low jitter frequency. It clearly shows that there is a 
series of asymmetric oscillations excited along the surface of liquid-vapor 
interface driven by asymmetric gravity jitter-dominated acceleration. 

Careful examination of the case of gravity jitter acceleration-dorainate 
liquid-vapor interface oscillations, gravity jitter acceleration, shown in 
Equations (5) and (6) which is also illustrated in Figure 3, indicates that there 
is a sinusoidal oscillation longitudinal to the direction aligned with spacecraft 
mass center to Earth center (parallel to unit vector r c shown in Figure 1). As 
we indicated in Equation (4), the angle xj> z varies with time. This phenomena 
shows that the gravity jitter acceleration exerted on the spacecraft is 
equivalent to time - dependent oscillatory forces which push the bubble in the 
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combined directions of down-and-up (see z-component of gravity jitter 
acceleration shown in Figure 3) and leftward- and- rightward (see x-and-y 
components of gravity jitter accelerations shown in Figure 3) as the bubble is 
rotating with respect to the spacecraft rotation axis. The down-and-up motion 
of the bubble is due to the fact that z-component gravity jitter acceleration 
decreases from the maximum downward direction to the minimum downward direction, 
during the time period between t - 0 to t - 300 s, which pushes the liquid 
downward and the bubble upward; while the gravity jitter acceleration changes to 
its downward direction and magnitude from the minimum upward to the maximum 
upward, and then back to the minimum upward, during the time period between t - 
300 s to t - 900 s, which pushes the liquid upward and the bubble downward; in 
the meanwhile, the gravity jitter acceleration changes to its downward direction 
and magnitude during the period from t - 900 s to 1200 s which repeats the half 
cycle story of the time period from t - 0 to t - 300 s. The lef tward-and- 
rightward oscillations of bubble are the results of positive and negative 
directions of gravity jitter accelerations in the x- and y-components as that 
illustrated in Figure 3. The positive direction gravity jitter acceleration 
implies that the liquid is pushed rightward and the bubble is driven leftward. 
The time variations of gravity jitter accelerations, shown in Figure 3, 
characterize the leftward and rightward oscillations of the bubbles shown in 
Figure 7. 

In addition to the 0.1 Hz low frequency gravity j itter -dominated 
acceleration, 1.0 Hz medium frequency and 10 Hz high frequency gravity jitter 
acceleration driven liquid-vapor interface oscillations have also been 
investigated. The characteristics of these three ranges of gravity jitter 
frequencies have been well-documented in the earlier studies 3 ' 6 . Because of page 
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limitation, oscillations driven by the time -dependent variations of the direction 
of gravity j itter-dominated acceleration with medium and high frequencies will 
not be shown in this paper. It is also worthwhile to mention that the natural 
frequency of a cryogenic helium container is much closer to the low frequency 
jitter acceleration than that of the medium and high frequencies. Results show 
that lower frequency gravity jitter-dominated acceleration contribute more in the 
driving of asymmetric profiles of the interface than the higher frequency gravity 
jitter-dominated acceleration. These results are agreeable with the earlier 
studies 3 ' 6 . 

IV. Characteristics of Slosh Wave Induced Fluctuations 
in Fluid System Moment and Angular Momentum 
Slosh wave induced fluctuations in the fluid system of the rotating dewar 
introduce time -dependent disturbances in moment and angular momentum of the 
spacecraft propulsion system. In this study, in addition to the existing angular 
velocity along the rolling axis, there are induced angular velocities along the 
yawing and pitching axes due to the sloshing dynamic modulated fluid motion 
inside the rotating container. These disturbances of angular velocities in the 
yawing and pitching axes readjust the angular velocity in the rolling axis. 

In order to accommodate the spacecraft dynamics of yawing, pitching and 
rolling, the cylindrical coordinates (shown in Figure 1) of the rotating 
container are transformed into cartesian coordinates based on (x, y, z) - (rco sfl, 
rsintf, z) with corresponding velocity components (V x , V y , V z ) - (ucos# - vsinfl, 
us infl + vcos# , w) . If the spacecraft is rotated with respect to the mass center 
(r c , $ ct z c ) in cylindrical coordinates, the location of the mass center in 
cartesian coordinates becomes (x c , y c , z c ) - (r c cos0 c , r c sin0 c , z c ) . The velocity 
of the mass center in cylindrical and cartesian coordinates which are defined in 
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the non- inertia frame is (u c , v c , w c ) and (v cx , v cy , v c 2 ) , respectively. 
Detailed mathematical derivations of induced fluid angular velocities, fluid 
angular momentum, and fluid moment modulated by sloshing dynamics driven by 
orbital accelerations are illustrated in our recent publications 16 * 21 and will not 
be repeated in this paper. 

(IV-A) Combined Gravity Gradient and 10~ 8 Background Jitter Accelerations 

Driven Angular Momentum Fluctuations 

As we mentioned earlier, the combined effects of gravity gradient and 
jitter accelerations are completely dominated by the gravity gradient 
acceleration when the 10‘ 8 ^ and less background gravity jitter accelerations 
are applied simultaneously with gravity gradient acceleration which drive fluid 
system moment and angular momentum fluctuations. 

Figure 8(A) shows the computed time variation of the fluctuations of 
angular momentum driven by gravity gradient-dominated acceleration. This figure 
shows the following results: (a) The values of angular momentum fluctuations are 
(AH X , AHy, AH Z ) - (2.765, 2.269, 5.094)10 5 g*cm 2 /s, it clearly indicates AH 2 > AH X 
> AHy. The maximum absolute values of angular moment are Max ( Kl, |h,|. | H j) 
- (1.57, 1.30, 4.166) 10 5 g*cm 2 /s. It also indicates |H Z | > |H X | > | Hy | . (b) 
The initial values of H x , Hy and H z start from zero in non- inertia frame while 
that of H z starts from a non-zero value in inertia frame. This is due to the 
fact that H z starts with a spinning speed of 0.1 rpm while there is no spinning 
velocity along x and y axes at the initial point in inertia frame. (c) 
Variations of oscillation frequency of H z are a reflection of spinning speed with 
0.1 rpm in which the peak values are at 300, 600, 900 and 1200 s and valley 
values are at 150, 750 and 1050 s. (d) Variation amplitudes of H x and Hy are 
mainly caused by gravity gradient acceleration and are smaller than that of H z 
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because there is no spinning motion applied on the x- and y- axes. 

Figures 9(A), 9(B) and 9(C) show variations of fluid moments due to 
sloshing dynamics driven by gravity gradient -dominated acceleration along the x, 
y and z axes, respectively. The values of fluid moment fluctuations are (AM*, 
AMy, AM 2 ) - (544, 526, 1085) dyne-cm. The maximum absolute values of fluid 
moment are Max (|mJ, |My| , |M,|) - (299, 285, 580) dyne cm. It shows AM* > AM, 
> AMy and |m,| > |M,| > | My | . Characteristics of the fluctuations of fluid 
moments driven by gravity gradient -dominated acceleration draw four point 
conclusions similar to those drawn for the fluctuations of angular momentum of 
fluid system shown earlier. 

Figure 10(A) shows time fluctuations of the locations of bubble mass 
centers of the rotating container due to slosh waves excited by gravity gradient- 
dominated acceleration. The values of bubble mass center fluctuations are (Ax c , 
Ay c , Az c ) - (4.54, 2.83, 2.32) cm. It shows Ax c > Ay c > Az c for bubble mass 
center fluctuations driven by gravity gradient-dominated acceleration. 

(IV) -B Combined Gravity Gradient and 10‘ 7 g Q Background Jitter Accelerations 

Driven Angular Momentum Fluctuations 

As we mentioned earlier, the combined effects of accelerations are equally 
weighted by both the gravity gradient and jitter accelerations when 10' 7 g Q 
background gravity jitter acceleration is jointly applied simultaneously with 
gravity gradient acceleration which drive fluid system moment and angular 
momentum fluctuations. 

Fluctuations of angular momentum, fluid moment and bubble mass center 
driven by combined acceleration for gravity jitter acceleration with frequencies 
of 0.1 , 1.0 and 10 Hz have been investigated. Because of page limitation, the 
results of fluctuations driven by combined acceleration for 0 . 1 Hz gravity jitter 
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frequency only will be presented in this study. Figure 8(B) shows the computed 
time variation of the fluctuations of angular momentum driven by equally weighted 
gravity gradient and jitter accelerations with frequency of 0.1 Hz. The values 
of angular momentum fluctuations are (AH X , AHy AH X ) - (3.266, 2.55, 10.308) 10 5 
g cm 2 /s and the maximum absolute values of angular moment are Max (|H X |, j Hy | , 
|H X |) - (1.666, 1.423, 9.472) 10 5 g cm 2 /s. It shows AH X > AHj > AHy and | H z | > 
|H X | > I Hy | • The results are very much similar to the four points conclusions 
driven for the angular momentum fluctuations driven by gravity gradient -dominated 
acceleration. It also shows that equally weighted gravity gradient and jitter 
accelerations have excited greater magnitudes and fluctuations of angular 
momentum than those driven by gravity gradient -dominated acceleration. 

Figures 11(A), 11(B) and 11(C) show variations of fluid moments due to 
slosh waves excited by equally weighted gravity gradient and jitter accelerations 
with frequency of 0.1 Hz along the x, y and z axes, respectively. The values of 
fluid moment fluctuations are (AM*, AMy, AM 2 ) - (838, 800, 1797) dyne cm and the 
maximum absolute values of fluid moment are Max (|m*|, | My | , |M 2 |) - (485, 468, 
965) dyne cm. It shows AM Z > AM X > AMy and |M 2 | > |M*| > | My | . Characteristics 
of the fluctuations of fluid moments driven by equally weighted gravity gradient 
and jitter accelerations draw four point conclusions similar to those drawn for 
the fluctuations of angular momentum of fluid systems driven by gravity gradient- 
dominated acceleration shown earlier. It also shows that equally weighted 
combined gravity gradient and jitter acceleration have excited greater magnitudes 
and fluctuations of fluid moment than that driven by gravity gradient- dominated 
acceleration . 

Figure 10(B) shows time fluctuations of the locations of bubble mass 
centers of the rotating container due to slosh waves excited by equally weighted 
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gravity gradient and jitter accelerations. The values of bubble mass center 
fluctuations are (Ax c , Ay c , Az c ) - (7.63, 6.27, 11.48) cm. A comparison of 
Figures 10(A) and 10(B) for bubble mass center fluctuations driven between 
gravity gradient-dominated acceleration and equally weighted gravity gradient and 
jitter accelerations provide the following conclusions: (a) Torsional moment-and 
twisting- force equivalent gravity gradient -dominated acceleration exerted on the 
rotating dewar container produces smaller values and fluctuations of angular 
momentum and fluid moment than that driven by equally weighted gravity gradient 
and jitter accelerations. (b) Dynamics of bubble (liquid-vapor interface) 
driven by torsional moment-and- twisting force-equivalent gravity gradient- 
dominated acceleration produces bubble mass center fluctuations of Ax c > Ay c > 
Az c while up and down oscillations of bubble driven by equally weighted gravity 
gradient and jitter accelerations induces bubble center fluctuations of Az c > Ax c 
> Ay c . (c) A comparison of the values of bubble mass center fluctuations shows 
that the fluctuations of bubble mass center induced by equally weighted gravity 
gradient and jitter accelerations are greater than those induced by gravity 
gradient-dominated acceleration. (d) Fluctuations of both x c and y c start from 
zero while that of the z c starts from a non- zero value at the middle point of the 
height of partially filled liquid container. 

(IV-C) Combine Gravity Gradient and 1CT 6 g Q Background Jitter Accelerations 

Driven Angular Momentum Fluctuations 

As we mentioned earlier, the combined effects of gravity gradient and 
jitter accelerations are completely dominated by the gravity jitter acceleration 
when the 10' 6 g Q and higher background jitter acceleration are applied 
simultaneously with gravity gradient acceleration which drive fluid system moment 
and angular momentum fluctuations. 
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Figure 8(C) shows the computed time variation of the fluctuations of 
angular momentum driven by gravity jitter-dominated acceleration with frequency 
of 0.1 Hz. The values of angular momentum fluctuations are (AH*, AHy AH Z ) - 
(16.013, 6.753, 29.164) 10 5 g cm 2 /s and the maximum absolute values of angular 
moment are Max (|H*|, | Hy | , |H Z |) - (8.667, 3.986, 29.165) 10 s gcm 2 /s. It shows 
AH Z > AH* > AHy and |H Z | > |H*| > | Hy | . A comparison of Figures 8(A), 8(B) and 
8(C) for angular momentum fluctuations driven by gravity gradient-dominated 
acceleration and gravity jitter-dominated acceleration shows that gravity 
gradient-dominated acceleration have excited smaller magnitudes and fluctuations 
of angular momentum than that driven by gravity jitter-dominated acceleration. 

Figures 12(A), 12(B) and 12(C) show variations of fluid moments due to 
slosh waves excited by gravity jitter-dominated acceleration with frequency of 
0.1 Hz along the x, y and z axes, respectively. The values of fluid moment 
fluctuations are (AM^, AMy, AM 2 ) - (2373, 2346, 5423) dyne’cm and the maximum 
absolute values of fluid moment are Max (|M X |, | My | , |m z |) - (1282, 1235, 3043) 
dyne cm. It shows AM Z > AM* > AMy and |M 2 | > Im^I > |My| . Characteristics of 
the fluctuations of fluid moments driven by gravity jitter-dominated acceleration 
draw four point conclusions similar to those drawn for the fluctuations of 
angular momentum of fluid systems driven by gravity gradient -dominated 
acceleration shown earlier. It also shows that gravity gradient-dominated 
acceleration have excited smaller magnitudes and fluctuations of fluid moment 
than that driven by gravity jitter-dominated acceleration. 

Figure 10(C) shows time fluctuations of the locations of bubble mass 
centers of the rotating container due to slosh waves excited by asymmetric 
gravity jitter-dominated acceleration. The values of bubble mass center 
fluctuations are (Ax c , Ay c , Az c ) - (17.28, 12.2, 33.83) cm. A comparison of 
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Figures 10(A) , 10(B) and 10(C) for bubble mass center fluctuations driven between 
gravity gradient-dominated and gravity jitter-dominated accelerations provides 
the following conclusions: (a) Torsional moment-and twisting force-equivalent 
gravity gradient-dominated acceleration exerted on the rotating dewar container 
produces smaller values of magnitudes and fluctuations of angular momentum and 
fluid moment than that driven by gravity jitter-dominated acceleration. (b) 
Dynamics of bubble (liquid- vapor interface) driven by torsional moment-and- 
twisting force-equivalent gravity gradient acceleration produces bubble mass 
center fluctuations of Ax c > Ay c > A z c while up and down oscillations of bubble 
driven by gravity jitter-dominated acceleration induces bubble center 
fluctuations of Az c > Ax c > Ay c . (c) A comparison of the values of bubble mass 
center fluctuations show that fluctuations of bubble mass center induced by 
gravity jitter-dominated acceleration is greater than that induced by gravity 
gradient- dominated acceleration. (d) Fluctuations of both x c and y c start from 
zero while that of the z c starts from a non-zero value at the middle point of the 
height of partially filled liquid container. (e) Fluctuation magnitudes of x c „ 
y c and z c associated with higher frequency jitter are lower than those of the 
frequencies associated with lower jitter frequency for bubble mass center 
fluctuations driven by gravity j it ter - dominated acceleration. 

V. Discussion and Conclusion 

The dynamical behavior of cryogenic liquid systems in particular, capillary 
effect governed liquid-vapor interface disturbances affected by the asymmetric 
combined gravity gradient and jitter accelerations at three different ranges of 
background gravity, with an example applicable to the GP-B Spacecraft dewar have 
been carried out by numerically computing the three-dimensional Navier- Stokes 
equations subject to the initial and boundary conditions. As gravity gradient 
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acceleration is pre -determined if the orbit of the spacecraft is fixed, the 
ranges of background gravity jitter acceleration will be different with various 
ranges of environmental drag due to the sources of residual acceleration, solar 
wind environment during solar minimum and maximum cycles, etc. Gravity jitter 
accelerations with three ranges of background gravity at 1CT 6 , 1(T 7 and 
10~ 8 go are chosen in this study. 

Induced angular velocities along the rolling, yawing and pitching axes 
caused by slosh wave excitation driven by gravity gradient -dominated and gravity 
jitter-dominated accelerations have been investigated. Results show that there 
is a series of large amplitude fluctuation of angular momentums and fluid moments 
along the rolling, yawing and pitching axes. In particular, slosh waves excited 
by the asymmetry gravity gradient-dominated and gravity jitter-dominated 
accelerations contribute greatly to the large amplitude fluctuations of angular 
momentum and fluid moment along the rolling axis of the spacecraft. The study 
of bubble mass center fluctuations driven by gravity gradient- dominated and 
gravity j itter-dominated accelerations show that gravity gradient-dominated 
acceleration drives greater bubble mass center fluctuations in the pitching axis 
while gravity j itter-dominated acceleration induces greater bubble mass 
fluctuations in the rolling axis. 

The time -dependent variations in the direction of gravity gradient and 
gravity jitter-dominated accelerations imposed on the spacecraft will change and 
modify the asymmetry distribution of liquid-vapor in the rotating dewar 
container. This means that fluid moment, angular momentum and bubble mass center 
will be deeply disturbed by the presence of asymmetric fluctuations in liquid- 
vapor interface. In this study, the differences in gravity gradient- dominated 
acceleration and gravity j itter- dominated acceleration induced slosh waves and 
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how these slosh wave excitation affected the fluctuations in angular momentum, 
fluid moment and bubble mass center disturbances have been studied. A better 
understanding of sloshing dynamics modulated fluctuations in angular momentum, 
fluid moment and bubble mass center distributions affected by different levels 
of orbital accelerations are extremely important for the development of 
sophisticated attitude and guidance control techniques to assure the proper 
operation of very complicated scientific spacecraft such as GP-B 22 . 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Coordinate system for the computation of gravity gradient 
acceleration, 0 S is the location of geometric center of cylinder. 
Direction of r c points toward Earth center. (x* , y* , z’) denotes 
absolute coordinate system, while (x, y, z) illustrates non- inertia 
frame spacecraft bound coordinate system. 

Time variation of gravity gradient acceleration acting on fluid mass 
located at (r, 8 , z) - (40 cm, ir/4 p 10 cm) for turn-around period of 
1200 s with rotating speed of 0.1 rpra. (A) Along x-direction; (B) 
Along y-direction; (C) Along z-direction. 

Time variation of gravity jitter acceleration acting on the fluid 
system under background gravity of 10" 8 , 10” 7 and 10' 6 g ol rotating 
speed of 0.1 rpm, turn-around period of 1200 s and gravity jitter 
frequency of 0.1 Hz. (A) Along x-direction; (B) Along y-direction; 
(C) Along z-direction. 

Initial profiles of liquid-vapor interface for Gravity Probe-B module 
of rotating dewar under background gravity of 10' 7 g OJ rotating speed 
of 0.1 rpm and direction of background gravity at — 0°. (A) In 

r-z plane at 6 - 0° and 180°, (B) In r -6 plane at 8 - 90° and 270°, 
(C) In r -6 plane at z - 108 cm, and (D) Three-dimensional liquid- 
vapor interface profile. 

Time sequence evolution of three-dimensional liquid-vapor interface 
oscillations for rotating dewar driven by combined gravity gradient 
and jitter accelerations with background gravity of 10~ 8 g 0 . 

Rotating speed of dewar is 0.1 rpm, and the period of gravity 
direction turn-around time is 1200 s. 
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Figure 6 Time sequence evolution of three-dimensional liquid-vapor interface 
oscillations for rotating dewar driven by combined gravity gradient 
and jitter accelerations of 0.1 Hz frequency, background gravity of 
10‘ 7 So with the period of gravity direction turn-around time - 
1200 s, and rotating speed of 0.1 rpm. 

Figure 7 Time sequence evolution of three-dimensional liquid-vapor interface 
oscillations for rotating dewar driven by combined gravity gradient 
and jitter accelerations of 0.1 Hz frequency, background gravity of 
10~ 6 g Q with the period of gravity direction turn-around time - 
1200 s, and rotating speed of 0.1 rpm. 

Figure 8 Time sequence of angular momentum fluctuations (H*, Hy, H z ) caused 
by the slosh wave excitation in rotating dewar with speed of 0.1 rpm 
and gravity turn-around time of 1200 s. (A) For angular momentum 
fluctuations driven by combined gravity gradient and jitter 
accelerations with background gravity of 10" 8 g G . (B) For angular 
momentum fluctuations driven by combined gravity gradient and jitter 
accelerations of 0.1 Hz gravity jitter frequency and background 
gravity of 10“ 7 g Q . (C) For angular momentum fluctuations driven 
by combined gravity gradient and jitter accelerations of 0.1 Hz 
frequency and background gravity of 10’ 6 g Q . 

Figure 9 Time sequences of fluid moment fluctuations caused by slosh wave 
excitation driven by combined gravity gradient and jitter 
accelerations of 0.1 Hz frequency and background gravity of 10" 8 g Q 
in rotation speed of 0.1 rpm and gravity turn-around time of 1200 s. 
(A) For moment along x-direction; (B) For moment along y-direction; 
and (C) For moment along z-direction. 
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Figure 10 Time sequences of bubble mass center fluctuations caused by the slosh 
wave excitation under gravity turn-around time of 1200 s and rotation 
speed of 0.1 rpm. (A) For bubble mass center fluctuations driven 
by combined gravity gradient and jitter accelerations with background 
gravity of 10' 8 g^ (B) For bubble mass center fluctuations driven 
by combined gravity gradient and jitter accelerations of 0.1 Hz 
frequency and background gravity of 10" 7 g^,; and (C) For bubble mass 
center combined gravity gradient and jitter accelerations of 0.1 Hz 
gravity jitter frequency and background gravity of 10' 6 g^ . 

Figure 11 Time sequences of fluid moment fluctuations caused by slosh wave 
excitation driven by combined gravity gradient and jitter 

accelerations of 0.1 Hz frequency with background gravity of 10‘ 7 g^ 
and rotation speed of 0.1 rpm under gravity turn-around time of 
1200 s. (A) For moment along x-direction; (B) For moment along y- 
direction; and (C) For moment along z-direction. 

Figure 12 Time sequence of fluid moment fluctuations caused by slosh wave 
excitation driven by combined gravity gradient and jitter 

accelerations of 0.1 Hz with background gravity of 10' 6 g or and 
rotation speed of 0.1 rpm under gravity turn-around time of 1200 s. 
(A) For moment along x-direction; (B) For moment along y-direction; 
and (C) For moment along z-direction. 
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INITIAL PROFILES OF LIQUID HELIUM AND VAPOR INTERFACE 
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(D) Three-Dimensional Profiles 
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